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ABSTRACT 


This study presents the derivation of a diffuser design parameter, 
which takes into account losses due to the geometry of the diffuser. 

A criterion for design of annular diffusers, using the design 
parameter to predict losses, is developed from reference data and 
model tests. 

A curved wall annular diffuser was designed and tested to verify 


the predictions of the design criterion for an arbitrary diffuser. 











1. 


ecT.. 


TV . 


Nis 


VI. 


Wo. 


TABLE OF CONTENTS 


INTRODUCTION 
DEVELOPMENT OF THE SHAPE PARAMETER OMEGA (Q) 

A. ENTROPY CHANGE OF A ONE-DIMENSIONAL, FRICTIONAL FLOW 
B. ENTROPY CHANGE FROM THERMODYNAMIC CONSIDERATIONS 

C. DEFINITION OF OMEGA (Q) 

D. RELATIONS BETWEEN Q AND DIFFUSER PERFORMANCE 

TEST INSTALIATION AND INSTRUMENTATION 

A, AIR DELIVERY SYSTEM 

B. DESCRIPTION OF DIFFUSER MODELS AND ASSEMBLY 

C. INSTRUMENTATION 

COMPUTATIONS AND EXPERIMENTAL PROCEDURES 

A. PRESSURE 

B. TEMPERATURE 

C. FLOW RATE 

D. CHECK OF NOZZLE FLOW RATE WITH SURVEY DATA 

E. REYNOLDS NUMBER 

ANALYSIS OF TEST DATA FROM REFERENCE 2 

DISCUSSION OF TEST DATA 


CONCLUSIONS AND RECOMMENDATIONS 


APPENDIX A - Determination of Diffuser Efficiency . 


APPENDIX B - Design of Diffuser Model (2) 


APPENDIX C - Computer Program Omega 


APPENDIX D - Reduction of Data from Sovran and Klomp, Reference 2 


Cl - Assembly of Input Data 


C2 - Output 


1) 


16 


i 
21 
Z2 
26 
26 
28 
33 
38 
Sie 
38 
39 
40 
44 
45 
48 
60 
61 
63 
67 
67 
68 


The. 











BIBLIOGRAPHY 


FORM DD 1473 


80 


83 








Table 


Ee 


ino 


Cl 


C2 


D1 


D2 


LIST OF TABLES 


Diffuser descriptions 
Model (1) and Model (2) test data 


Meridional static pressure values for 
Model (1) and Model (2) 


Listing of Program Omega 
Sample output of Program Omega 
Listing of reference data reduction program 


Sovran and Klomp, GMR-511, data 


Page 
35 
il 


58 


70 
74 
76 


77 








Figure 


10 
Wal 
12 
13 


14 


15 


16 
17 


18 


19 


LIST OF ILLUSTRATIONS 


Example of diffusers which previously had 
the same description 


One-dimensional flow in a diverging channel 


T-s diagram for an adiabatic, polytropic 
compression in a stationary duct 


Schematic of the air delivery system 
Transition ducting, looking downstream 
Support assembly, including support struts 


Diffuser support assembly and converging 
outer wall (with Model (1) attached) 


Diffuser support assembly with Model (1) 
cantilevered center-body installed 


Center-body of Model (2) 


Outer wall of Model (2) with support assembly 
Schematic of diffuser assembly 

Schematic of diffuser Model (1) 

Overall view of the experimental apparatus 
Diffuser inlet plane, velocity survey in 
radial direction (measured by hot-wire 
anemometer ) 

Diffuser inlet plane, velocity survey in 
radial direction (measured by five-hole 
pressure probe) 

er vs. Q for Reference 2 data 


ce vs. 9, Model (1) test data 


Comparison of Reference 2 data with 
Model (1) and Model (2) test data 


Diffuser polytropic efficiency vs. 
Reynolds number 


Page 


14 


17 


20 


27, 
29 
29 


30 
30 


31 
ol 
32 
34 
36 


41 


42 


46 
49 


50 


54 








20 


Zi 


Ze 


ZS 


24 


Dimensionless pressure distribution along 
outer wall of Model (1) 


Dimensionless pressure distribution along 
the inner and outer walls of Model (2) 


Design of Model (2) 
Pichine Of circular arcs 


Input field specifications for Program Omega 


55 


56 


64 
65 


69 








Symbol 
A 


C 


rl 


~~ | atl 


™y 


mV 


Ap 


List OF sox tBOLS 


Definition 

Area 

Wetted perimeter 

Coefficient of pressure recovery 

Ideal coefficient of pressure recovery 
Local coefficient of friction 

Average local coefficient of friction 


Specific heat of air at constant 
pressure 


Distance from outer wall 
Unit conversion factor 


Texas Insirument pressure gauge 
Conversion Factor 


Axial length 


Meridional length along mean 
streamline 


Meridional wall length 
Meridional length along wall 
Mach number 

Mass flow rate 

Referred mass flow rate 
Temperature reading 
Polytropic exponent 

Total pressure, absolute 
Static pressure, absolute 


Pressure difference between two 
points in the flow 


Unies 
in 


in 


in 

2 

ft-lb /lb-sec 
m 


in H, 0/count 


in 


in 


in 


in 
slugs/sec 


millivolts 
haya 
Wits 


sy is 








= 


Subscripts 


1 


lr 


Dynamic pressure (P. - p) 
Gas constant 

Reynolds number 

Entropy 

Total temperature 
Static temperature 
Temperature 

Temperature 

Average velocity 

Flow rate 

Equivalent flow rate 
Factor (P, - Pgd/Pa 
Ratio of specific heats 


Density coefficient 


Dimensionless static pressure 


Polytropic efficiency 
Dimensionless length 
Viscosity 

Density 

Wall shear stress 
Flow function 


Shape parameter 


Station at diffuser inlet plane 


Referred to station one 


Station at diffuser exit plane 


Atmospheric 


10 


ina 


ftv-co bb /lbo Scr 
m 


ft/sec 


lb /sec 
m 


lb Be OvIETa 
m 

3 

slugs/ft 


ib/£t? 








TI 


Downstream from nozzle 
Hub or inner wall 
Isentropic 

Poly Grop 1c 

Tip or outer wall 


Values measured on Ghee texas Instrument 
Fused Quartz Pressure Gage 


Total 
Upstream from nozzle 


At or along wall 


11 








AC KNOW LEDGEMENT 


The author wishes to acknowledge the guidance given by Professor 
M. H. Vavra of the Department of Aeronautics of the Naval Postgraduate 
School, and assistance of Mr. James E. Hammer and the technical staff 
of the Turbopropulsion Laboratory in the preparation of the experi- 
mental portion of this thesis. 

Special thanks go to Mr. Gordon Gulbranson, non-metal model maker 
of the Department of Aeronautics for the beautiful diffuser models 


which he fabricated. 


iz 








I. INTRODUCTION 


Im the design of annular diffusers, such as those found in jet 
engines, and turbomachines, there has been a problem defining param- 
eters which indicate sufficiently the optimum design range for varying 
geometries. Many schemes have been devised, such as the so-called 
"equivalent conical angle method", to relate a diffuser to a known 
shape. The difficulty in these approaches has been that each person 
using the method developed a slightly different approach which resulted 
in a wariety of descriptions for the same diffuser. Gleason [1] gives 
a summary of such methods and points out the difficulty in correlating 
data from different sources. 

Some apparent problems arise when describing a diffuser by inlet 
area, exit area, length, difference in inlet radii, equivalent cone 
angle, or combinations of these parameters developed by empirical means. 
Although these terms give some intuitive feel for the way a diffuser 
would perform, two diffusers which differ widely in performance could 
have the same definition. 

The simple example of Figure 1, where the two diffusers have the 
same inlet area, exit area, length and possibly the same divergence 
angle, makes it immediately obvious that the two diffusers would not 
have the same losses. The velocity and pressure distributions through 
the diffuser channel have a direct effect on the boundary layer 
development, which in turn controls the frictional losses. The 
widely differing geometries of the two above-mentioned diffusers 


would cause large differences in boundary layer growth. Therefore, 
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FIGURE 1 


Example of diffusers which previously had the 
same description. 
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the performance of the two would differ and a similar description 
would not be appropriate. 

Sovran and Klomp [2] have empirically developed a fairly successful 
set of parameters for straight wall annular diffusers. The geometric 
area ratio was taken as the primary factor affecting diffusion, and a 
non-dimensionalized length when combined with the area ratio, gave an 
overall criterion for the pressure-gradient affecting the boundary 
layer development. Their description of annular diffusers is at 
present limited to designs with straight walls. 

Vavra [3 and 4} presented the idea of a shape factor or parameter 
which includes the effects of diffusion and boundary layer influence 
for the evaluation and preliminary design of diffusers of arbitrary 
shapes. 

The objectives of this study were to present the derivation of the 
shape parameter, 9, to substantiate its usefulness from reference data, 
and to evaluate the validity of the results by construction and investi- 
gation of a family of straight wall diffusers such as those in Ref. 2. 
To show that the design shape parameter may be applied to diffusers 
other than those with straight walls, a curved wall diffuser such as 
might be found at the exit of the high pressure compressor of a jet 


engine was designed for an optimum () and evaluated by tests. 
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II. DEVELOPMENT OF THE SHAPE PARAMETER OMEGA ((Q) 


The function of a diffuser is to convert a portion of the kinetic 
energy of the free stream into static pressure by a compression process 
similar to that in a mechanical compressor. The shape parameter, 

Omega (), to be used for diffuser design, can be evolved by examining 
the differential entropy changes of such a polytropic compression 


process. 


A. ENTROPY CHANGE OF A ONE-DIMENSIONAL, FRICTIONAL FLOW 

Am expression for the differential entropy change of a polytropic 
compression may be found by investigating a one-dimensional flow in a 
diverging channel, such as shown in Fig. 2. A momentum balance on the 
differential volume of Fig. 2 gives 


m (V f dV) - mV = pA - (p # dp) (A # dA) f (p # dp/2) dA - TCdL 


which reduces to oh) 
f 
m dV = - Adp - TCdL (2) 
or 
= 
dp = ii dv A dL (3) 


where the values of A, C, p, dL, V, and T are as defined in Fig. 2. 
Introducing the equation of continuity, 
m = oVA (4) 


into equation (3) results in 


dp = - pvay - 1504 (5) 
OF 

82 52 a0) =Lo au 6 

é wy2) -2§ (6) 
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FIGURE 2 


One-dimensional flow in a diverging channel. 
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with the First Law of Thermodynamics 


dq du # pdv = dh - vdp = dh - - dp = Tds (7) 


or 


“~ dh - Tds (8) 
Equating (6) and (7) gives 


dh # d(v7/2) - Tds = dL (9) 


DjA 
PIA 


The first two terms of (9) combine to the derivative of the total 
enthalpy, H, and if the flow is restricted to an adiabatic case, with 
no energy in any form transferred between the fluid and its surroundings, 
the total enthalpy is constant and its derivative is equal to Zero. 
Therefore, 
dH =d (h # V“/2) = 0 


and equation (8) simplifies to 


Tds = 


D [A 
10 


dL (10) 


Introducing the concept of the local coefficient of skin friction 


in @ one-dimensional flow, 


= fe Sa la 
1/2pV 


equation (10) becomes 


2 


“¢ 
Tds Tar wa alll (12) 


C 
A 
Introducing equation (4) and the equation of state for a perfect gas 


in (12) results in an expression for the differential entropy change, 


Oy 2 oe 
T 2 ¢ 





nel le 


C 
—z dL (13) 
A 
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B. ENTROPY CHANGE FROM THERMODYNAMIC CONSIDERATIONS 
The differential change in entropy can also be found for an 
adiabatic, polytropic compression process as shown in the T-s diagram 


of Fig. 3. The definition of polytropic efficiency is, 


cee 

"Np “aT oe 
The isentropic relation for a differential compression from p to 
p + dp is 

y-1 
ih Ge ralth y 
= eS: (2+ se (15) 
T p 


By expanding (dp/p) by a binomial series, and assuming dp to be small 


compared to p, 
dT mee (16) 


is Y 


Combining (14) and (16) gives 
qp - L 2! op (17) 
lee) - oP 


Solving equation (7) for 


ds Sue Go oe a 
and introducing 

“p ” Re yl 
and the perfect gas law 

p/o = R, 7 


results in an expression for 


a x af _ dp] 
wer (2 -s (18) 


subscetuting (17) into (18) gives 


ds =R C125 ee 1| (19) 
& P Np 


ED 











FIGURE 3 


T-s diagram for an adiabatic, polytropic compression in 
a stationary duct. 
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C. DEFINITION OF OMEGA (() 
Equating the two values of differential entropy change from 


equations (13) and (19) results in 
Cee R R 
sf RS Bay = -B[L - 1] ap (20) 
2 pe 3 P Ptoi 


Introducing the concept of referred flow rate, a dimensionless quantity, 


R, T 
m = 
tr Py] 


- 





m (21) 


and substituting (21) into (20) and rearranging gives 


e A 2 
st ty x) f ap=[+-i1]2a(2) (22) 
1 Np 1 1 
For a polytropic process 
n 
p/p = constant 


where n is the polytropic exponent. 


Rewriting (22) in the form 


3 1/n 
Se . RINE 1 
2G) Cae[L-]@) 2 as 
1 Np Py Py 
gives two independent differentials which can be integrated. 


Defining the left side of (23) as dx, and integrating, 
3 


Cc 
Eta Lae & 
fa? Gl a. al (24) 
l 


where L is the meridional length along the mean streamline. 


x, = 


oes 


Assuming an average value of coefficient of friction, Que for the 
whole surface and that my, is a constant for a given flow rate, (24) may 


be rewritten as 
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ee | 


ia SN 
Kaa ae (ty ¢ - al (25) 


" 


The imtegral of (25) is a parameter which contains the flow area ratio 


the diffuser length, and the wetted surface, and will be called the 


"diffuser shape parameter", defined by 
iL 3 
Q = f jG) eerie (26) 
0 


D. RELATIONS BETWEEN (1) AND DIFFUSER PERFORMANCE 


Defining the right-side of (23) as dX 


23 
Po/ Py ; 1/n 
X, = f [> - 1| p_) a (27) 
ne Py Py 


If Np and n are assumed constant, 


CUA), 
2 = begs oe Cine m) , 1D IG aa 1] Co, 


PS 
{ 


For 
Po = Py + Ap 


the binominal expansion 


@ l/n + 1 p, + A 


1/n + 1 l/n + 1 
Ge vam, Ger, 


1 +4 1) SP 4 . 
Py 


Assuming Ap is much less than Py and omitting higher order terms from 


the binomial expansion, equation (28) becomes 


x, =(+-1|2 (29) 
Ny Py 


Ze 








Equating xy and X 


Pig 


Q | 


ie See 8 Be 
Py J IE 


(30) 
With the equation of continuity, equation (21) becomes 
JR T 
mee ie ann er 
JEseJl | 
Substituting my, into equation (30) 
|= - 1| —P. = 9 (32) 
Np 1/20,V Vy 
For a constant total temperature, equation (18) becomes 
dP 
t 
ds = - ae P, 
Thus equation (19) becomes 
dP 
t ur P 


Integrating (33) 


Ge) L-]n@ 


(34) 

Using the assumption that P_ = Ee + AP. and Po = Py + Ap, where AP 
2 1 

and Ap are small compared to P. and p, respectively, and with the 


truncated power series expansion of In (1 + x) = x, (34) becomes 


AP, 
ty "p Py 


or 


a ey, (35) 
Paap 2M S 


1a) 








Za 
Since for an incompressible flow, M_, is small compared to unity, 
L 
| - 1| Ap = AP (36) 
Ys ic 


Substituting (36) into (32) gives 


7 AP, 
6. ee (37) 


7 
W2 P, Vy 


In Ref. 2, the Pressure Recovery Coefficient is defined as the ratio 
of the actual pressure rise and the maximum attainable one for the 


particula: flow rate, or 


PP 
¢ eel a 
P Wy pv, 


Defining the Ideal Pressure Recovery Coefficient as that which 
would result for an isentropic compression producing the theoretical 


pressure py of Fig. 3. 
th ry 
= J eh : 
iy: G9) 
P 1/2 pV, 


For isentropic flow AP. = 0, and assuming density to be constant and 
introducing the equation of continuity, 


0/2 v, - 0/2 V5" AY A 
GS, = =] - —) (40) 
pri 2 


p/2 Vv, 





From equation (40), 


2 


2 
p/2 V, 


Then 


2 2 
=p, t 0/2 V) - Py - p/2V, f1-c ] 


ae) 
1 

rd 
| 


pri 


Or 








t t eee 
p/2 Vv, p/2 V, P 


Using the definition of wa (38), and equation (37) there is 


NC ¢ = ei ~ oe (42) 
or 

= C oe C 

By = ea (43) - 


Hence , after having previously determined (), the average frictional 


loss coefficient, ¢,, can be calculated from the measured quantities 


Py> qy> Ais Po » and A, 


72) 








III. TEST INSTALIATION AND INSTRUMENTATION 


The apparatus used in this investigation was located in Building 
230 at the Turbopropulsion Laboratories of the Department of Aeronautics, 


Naval Postgraduate School, Monterey, California. 


m, AIR DELIVERY SYSTEM 

A schematic drawing of the installation is shown in Fig. 4. The 
air source was a 200 H.P.’, Carrier centrifugal compressor with a maxi- 
mum flow rate of about four pounds mass of air per second at a pressure 
ratio of about two. Gross throttling was accomplished by a main by-pass 
valve ahead of the air coolers. Fine adjustments of the flow rate were 
made by bleeding air from plenum (1) through a secondary by-pass valve. 

The air from plenum (1) discharged through a flow straightener into 
a four inch pipe containing a flow measuring nozzle. The nozzle used 
was the same as that described by Kelly [5], which had a diameter of 
3.022 inches and a diameter ratio, B, of 0.734. 

The flow was diffused into plenum (2) by a ten degree conical 
diffuser extending downward from the inlet. The discharge of plenum 
(2) was into a rectangular wooden duct containing flow straighteners 
and a total temperature, iron-constantan thermocouple probe. 

The mean stream velocity in the duct at maximum flow rate was 
approximately 20 to 25 feet per second. Transition pieces shown in 
Fig. 5 were used to change the geometry of the flow channel from 
rectangular to circular and to accelerate the flow to the diffuser 


inlet velocity and pressure. 
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B. DESCRIPTION OF DIFFUSER MODELS AND ASSEMBLY 

The diffuser assemblies of the two models constructed and their 
supporting assembly are shown in the photographs of Figures 6 through 
10. The supporting assembly consisted of the duct transition pieces 
with a cylindrical portion, shown in Fig. 6, containing three equally 
spaced support struts for a cantilevered center-body, so that no 
obstructions or supports would be necessary in the diffuser flow 
channel. The flow was accelerated to the diffuser entrance in an 
annular space created by the nose of the center-body assembly and a 
converging outer wall. 

A straight section, where the annular area remained constant, was 
provided just downstream of the inlet to provide room for flow measure- 
ments and surveys, as shown in the schematic of Fig. ll. Smooth 
transitions to the diffuser contour were arranged as shown in Figures 
8 and 9. All computations of the diffuser shape parameter included the 
flow area from the inlet plane to the exit plane of the diffuser, 
inclusive of the cylindrical inlet section. 

The diffuser inlet plane was provided with three supports for 
probes, equally spaced around the outer periphery and bisecting the 
angles formed by the supporting struts of the center-body. These 
Supports could accommodate both hot-wire and pressure probes for 
radial surveys. One support was used for a Kiel probe for the measuring 
of total pressure at the inlet plane. Three static pressure taps were 
provided in the same plane, 120 degrees apart on the outer wall. Three 
sets of static taps were placed axially downstream from the inlet static 
taps, at two-inch intervals on the outer wall of the first model, and 


on both the inner and outer walls at one-inch intervals on the second 
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FIGURE 5 


Transition ducting, looking downstream. 





FIGURE 6 


Support assembly, including support struts. 
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FIGURE 7 


Diffuser support assembly and 
converging outer wall (with 
Model (1) attached). 





FIGURE 8 


Diffuser Support assembly with 
Model (1) cantilevered center-body 
installed. 
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FIGURE 9 


Center-body of Model (2) 





FIGURE 10 


Outer wall of Model (2) with Support Assembly. 
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model. Static pressure readings were found to be equal at all three 
stations around the circumference, and the three taps at each axial 
location were therefore connected together. 

Model (1) had straight walls as illustrated in Fig. 12, with a wall 
divergence angle of 15 degrees for both the inner and outer contours. 
The construction was such that a series of four diffusers with differ- 
ent shape parameters could be produced by simply cutting off the end of 
the model. Model (2) had a more complicated shape, and was designed 
as described in Appendix B to show the validity of the parameter for 
curved-wall diffusers. A summary of the geometrical characteristics 


of both models is given in Table I. 


C. INSTRUMENTATION 

Pressure readings for the flow nozzle and the diffuser inlet were 
made with a Texas Instrument Fused Quartz Pressure Gage, which is 
shown in ane overall view of the experimental apparatus of Fig. 13. 
A constant pressure source calibration of the Pressure Gage, to a 
Meriam Vernier Manometer filled with water, gave a constant value of 
0.800 inches of water per count, with a reading accuracy of + 0.001 
counts. This result was identical to that of Beck [6]. 

Pressure readings for five-hole probe measurements were made on a 
235 cm, eight tube water manometer bank, which can be read to an 
accuracy of + 0.1 cm. Static pressure measurements along the diffuser 
length were made on a 92 inch, 20 tube water manometer board with a 
reading accuracy of + 0.05 inch. All pressure taps were connected to 


their respective reading devices with plastic tubing. A standard 


mercury wall-barometer was used for atmospheric pressure readings. 
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Temperature readings from the thermocouples upstream of the flow 
nozzle, and in the wooden duct, were obtained with a Leeds and Northrup 
portable potentiometer with compensated cold-junction, giving a read- 
ing accuracy of + 0.01 mV. Adiabatic flow was assumed from the area 
of temperature measurement to that of pressure measurement. 

Hot-wire measurements were conducted with a two channel, constant 
temperature, hot-wire anomometer manufactured by Applied Science, 
Carmel, California. The wires were 0.00015 inch diameter tungsten with 
a 0.1 inch gap. Surveys were taken radially from the outer to the 


famner wall. 
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IV. EXPERIMENTAL PROCFRDURES AND COMPUTATIONS 


All reduction of data was carried out by the Fortran IV computer 


program ‘'Omega'’ described in Appendix C. 


A. PRESSURE 
Pressures obtained from the Texas Instrument Fused Quartz Pressure 
Gage were reduced to pounds per square inch by 


p = 0.03613 K py, (44) 


where K = 0.8, the conversion factor from counts to inches of water, 
Pry is the pressure read in counts on the Pressure Gage, and the 
constant value, 0.03613, is the conversion factor for inches of water 
to pounds per square inch from Ref. 7. 

Barometric pressure was converted from inches of mercury to pounds 


per square inch by equation (3-2) of Ref. 6. 


0.0001634% t 
c | 


1 + 0.0001818 oe (49) 


py = 0.4912 p.| 1 : 


where the coefficients of tO: the room temperature in degrees centi- 
grade, are the thermal expansion coefficients of mercury and brass. 
The factor 0.4912 is the conversion factor for inches of mercury at 


0°C to pounds per square inch from Ref. 7. 


B. TEMPERATURE 
Thermocouple readings, recorded in millivolts, were converted to 
degrees Fahrenheit by the temperature relationship, 


te = 32.144 # 35.77 mV - 0.4518 mV (46) 
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obtained from a least squares fit of data from Ref. 8. This relation 
holds for temperatures below 100 degrees Fahrenheit, with the cold- 
junction maintained at 32 degrees, or compensated to 32 degrees by the 


standard cell of the Leeds and Northrup measuring apparatus. 


C. FLOW RATE 

After setting the desired flow rate, the fluid temperatures were 
allowed to stabilize before data were taken. Upstream static pressure, 
Pi? differential pressure across the flow nozzle, (P-Py)> and the 


upstream total temperature, t are necessary to compute flow rate 


Fug 
as described in Appendix A of Ref. 5. 


p fe 
* u (a 
a ie 47) 
fee, Ve 
P - p 
mae) one ew (48) 
t Pa 


* , ; : ee 
where w is the equivalent flow rate, which has dimension in . 


and 


= 


The so-called expansion coefficient, ¢, is 


426 0.286 
6 = |2.433 nox) eae (49) 
.305(1-x) 
where 
eae hd 
PL (50) 
and 
te - “ty 
E. = 1 # 0.0015 Sn (51) 


is a correction factor that takes account of the thermal expansion of 


mae nozzle throat diameter. 
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D. CHECK OF NOZZLE FLOW RATE WITH SURVEY DATA 

Preliminary surveys at the inlet plane with a hot-wire probe 
indicated that the velocity profiles at the three stations surveyed 
were similar. Figure 14 is a least-squares fit of the data points from 
all three survey stations. The data were taken at varying flow rates 
at each station, and the velocities were found to vary within + 1 
per cent of the mean obtained with the least-square fit. Hence, the 
flow was very nearly axisymmetric, and the non-dimensional profiles 
remained unchanged at the different flow rates investigated. The 
boundary layer was very thin at the inlet. The actual thicknesses were 
0.025 inch on the inner wall and 0.75 inch on the outer wall according 
to the measurements. 

To verify the pressure and temperature measurements, a comparison 
was made of the flow obtained from the survey data at the inlet and 
that measured by the flow nozzle. A survey with a five-hole probe was 
made at the inlet plane to obtain the actual values of the velocities 
which are plotted in Fig. 15. Data for the five-hole probe were cor- 
rected by use of the calibration curves supplied by the manufacturer, 
United Sensor and Control Corporation. The velocity profile obtained 
was essentially the same as that of the hot-wire surveys. The hot-wire 
values were used for the boundary layer regions, where the immersion 
Pe eet ions for the pressure probe were unreliable. The yaw angles of 
Gie=r low were less than + eee and the flow had pitch angles between 
three and six degrees, with the velocities pointing outward from the 


center. The greatest velocities occurred at the outer wall. 
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The five-hole probe survey showed the total pressures and the 
velocities to be within two per cent of the mean values in the turbu- 
lent central region outside the boundary layers. An average Mach 
number of the flow in this central region was determined from the 


isentropic relationship 


y-1 1/2 
ee J 
{FAlG - uh 
re Sy —, L (52) 
The average static temperature at the inlet was determined by 


assuming adiabatic flow from the upstream total temperature probe to 


the inlet plane, which gave the isentropic relationship 


T 
ot tigers ein (53) 
T rol 

1 

Hence, 
T 
tl 

ee (54) 

Saveset: 


Static pressure was found to be nearly uniform across the inlet 
plane, so that density was almost constant at this station. From the 


perfect gas relationship, 


p, (144) 


o, =————— 
1 S, g. qT, 


(55) 


with the units of slugs per cubic foot. 


From the survey, the velocities were known as a function of radius, 


V(r). Thus, a flow rate was calculated at the inlet plane by numerical 


integration of 


Ry 
Wy =2T Py f Vitis eo ae (56) 


os 
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i= = 
—_—— = 
— a 





assuming axisymmetric flow and no blockage from support strut wakes. 
This flow rate was found to be three to four per cent higher than the 
flow rate obtained from the flow nozzle measurements. This comparison 


was considered satisfactory for the range of measuring capability. 


E. REYNOLDS NUMBER 
The characteristic length for the Reynolds number (Re) was taken 


as the hydraulic radius of the inlet annulus, AR,,as in Ref. 2, or 


i ts 
2 (flow area), 


aT i (wetted perimeter), SD, 
Therefore, 
pe Vea 
| ee | i 
Ree Ss 58 
1 (58) 
and, substituting 
k 2 2 
© =p, V,4,=0,V,7 R-R) (59) 
I 1 


the Reynolds number (Re) can be expressed in dimensionally correct units 


as 


___w (12) 
Re = oun pF Ry) soe) 


f lb Layey SEES 


where viscosity, UW, has a representative value of 3.9 x 10) £ 


for the operating range of 85 to 95 degrees Fahrenheit. The flow rate, 


w, is in Ib /sec, and the radii are expressed in inches. 
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V. ANALYSIS OF TEST DATA FROM REFERENCE 2 


Sovran and Klomp [2] made a study of more than one hundred annular 
diffuser configurations, and their published data was of a form suit- 
able to be analyzed for the present purpose. Of the many other works 
published, none furnished a sufficiently wide cross-section of data 
which could be readily applied. Data from Ref. 2 were reduced by the 
Fortran IV program shown in Table Dl. The results are listed in Table 
D2 and plotted on Fig. 16 for a Reynolds number of approximately 6 x 10°. 
Each set of data points represents a family of straight wall diffusers 
having the same inlet area and wall divergence angles as shown in 
Fig. 16. The data points fall within a band on the plot of 7 versus 
Q, which is bounded by lines representing a dimensionless length 
description of the diffuser, L/AR. An investigation of the values of 
L/AR for existing diffuser designs indicates that the boundaries shown 
on Fig. 16 are in the range PW ceieai designs. The points falling 
below L/AR = 1.5 represent designs which become so short that diffusion 
is degraded, and those above L/AR = 13 become excessively long for 
practical use. It can be seen, that as () falls below approximately 


five, there is an abrupt rise in the value of Q¢,, which is repre - 


f? 
sentative of the entropy rise and thus of an increase of the losses in 
the diffuser. The trend of the curve for larger values of {Q was found 
by considering that, as a diffuser with small divergence angles becomes 
long, its characteristics approach t.iat of a pipe. Reference 9 gives 


the value of the coefficient of friction for an annular pipe as 


e. = 0.01 at Re = 6 x 10°. A pipe with L/AR = 13 has a value of 
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FIGURE 16 


Q Gr. vs. Q for Reference 2 data. 





Q = 26 and ce = 0.26, which is a higher value of for than that 
occurring from diffuser tests at = 11.5 and L/AR = 13, suggesting 
an upward trend in the curve for high values of (). This trend 
indicates that a minimum of ace occurs between about ( = 5 and 
Q = 15, possibly at a value of ) = 10, for an L/AR = 13. 

Wayra [3] calculated values of Qe for two-dimensional subsonic 


diffusers, which lie within the same band as the data plotted on Fig. 


Hee 
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VI. DISCUSSION OF TEST DATA 


The test results obtained from diffuser Model (1) and Model (2) are 
given in Table II. The quantities Oe for Model (1) are plotted on 
Fig. 17 versus (). Tests were made at different flow rates to investi- 
gate the effect of Reynolds number on the ieee an of points in the 


-) oF vs. Q plot for the different configurations of Model (1). Figure 


17 shows that as Reynolds number increases the quantity Qe decreases, 


f 


and vice versa. This effect could be anticipated from the trend 3: On 


against Re for flows in annular pipes, as shown in Ref. 9. 

The curve for L/AR = 6 and the curve for the family of diffusers 
with a = e. = ae. from Fig. 16, and the curves from Fig. 17 were 
transposed to Fig. 18. The shape of the curves from data of Model (1) 
approximates that of the similar family of diffusers from Fig. 16 for 
oF = e. = Se However, the values obtained from the tests of Model (1) 
are considerably lower than those from Ref. 2. This discrepancy may be 
explained by the fact that the models of Ref. 2 were constructed of wood 
and, therefore, had a higher wall roughness than the cast resin models 
used in this study. The value aE G, for a wood surface would be 
larger than that for a resin surface with lower wall roughness, similar 
to what can be seen from the Moody diagram for pipes of differing 
roughness. 


The values of Qc, for configuration four of Model (1) became 


f 


slightly negative. Values of (0 Si. from Ref. 2 data were also negative 


when Bis: became greater than a For an adiabatic process, the value 


Qe 


cannot be negative, since C is less than C __, by definition. An 
if pr pri 


investigation of the pressure meaSurements for configuration four of 
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TAREE TUL 


Model (1) and Model (2) test data 


Pi oZi 


iL 
2 9.689 
3 S20 tZ 


1 





roe ees CG S-  O.lC SO 


4 6.470 -0.0162 18209 lh QZ 
-0.0027 .8756 “> 
-0.0105 5S 7 1.034 
=-0'.0128 .8677 1.265 
=O toa .8589 1.494 





TABLE II, continued 
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Model (1) showed a fluctuation of the instrument readings of sufficient 
magnitude for C to be larger than C _,. 
pr pri 
It can be anticipated that the curves of Fig. 16 will be shifted 
vertically up or down with Reynolds number similar to the condition in 
Fig. 17. The values of ) for the point of sharply increasing © €. 
would not be affected by a change in Re, since the entire band of data 
points would be displaced in the vertical direction only. To avoid 


large values of Qc., a minimum value of Q = 5 should be used for most 


fee 
practical designs. 

To demonstrate the usefulness of the design parameter (9, for com- 
plicated shapes, Model (2) was designed as Meine in Appendix B for 
a value of Q = 5.7, and L/AR = 6.85. Values of Q = 6 and L/AR = 6 
were chosen for the original design point, to obtain an o Be = 0.1. 

By extrapolating from the curve for L/AR = 6 of Fig. 16 to the 
data of Fig. 17 to show the effect of Re and the Pe in the 
model surfaces, the heavy dashed curve of Fig. 18 was obtained. The 
data of Model (2), plotted on Fig. 18, verified the location of the 
L/AR curve for an arbitrary shaped diffuser. 

The efficiencies for both Model (1) and Model (2), plotted on Fig. 
19, remained essentially constant through the Re range of this study. 


Static pressure (p_) on the inner and outer walls of Model (1) and 
p P, 


Model (2) was non-dimensionalized by 


C2 sa . (61) 


and plotted against non-dimensionalized, meridional wall length, 
rN = 1/L, on Figures 20 and 21, where lis the meridional distance 


along the wall from the inlet to the pressure tap, and fig is the 
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meridional wall length. Values for C, A, and p,, are contained in 
Table III, for the outer wall of Model (1) and for both the inner and 
outer walls of Model (2). The static pressure became higher than the 
ambient pressure at about X} = 0.5 on the inner wall of Model (2), 
which gives a € less than that at the exit. The value of € reaches a 
minimum between 4 = 0.7 and A=0.9 and increases again toward € = 0.293 
at the exit plane. This condition may be due to a separation of the 
flow, and indicates the need for proper wall contouring in addition to 


satisfying the criterion of 9. 
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TABLE IIT 


Meridional, static pressure values for Model (1) and Model (2) 











CONTOUR (in. water) 
gage 

OUTER 

CONFIG. 1 

Model (1) 

OUTER 

CONFIG. 4 

Model (1) 


* } - The dimensionless length We, 


*kC - Dimensionless pressure 


58 





TABLE III, continued 


ee a eee 





CONTOUR STATION 
INNER 05107 0.787 
Model (2) 
0.558 
0.426 
Opis) 
0.280 
0.248 
On231 
0.241 
UG 
OUTER 1 0.0 -11.41 1.000 
Model (2) | 
2 Obs 74 1 -10.91 Ue io, 
3 0.241 -7.90 0.783 
4 0.366 -5.18 0.614 
5 0.484 -3.50 0.510 
6 0.605 -2 .40 0.442 
7 0.718 -1.58 Oneg! 
8 0.823 -1.15 0.364 
9 0.938 =O. 27 0.316 
10 F000 0.00 On293 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


This paper presents a design parameter (), for diffuser designs, 
which takes into account losses due to the geometry of the diffuser. 

It is shown from reference data and model tests that there is an 
abrupt increase in losses for diffusers designed for a value of (Q 
less than about five, and that this effect is independent of Reynolds 
number. 

The parameter ()} is defined for incompressible, adiabatic flows. It 
is recommended that further experiments be made to extend the data for 
incompressible flows to higher values of 0, and that an investigation 
of the compressible case be made. 

It is suggested that Model (2) be improved by recontouring the 
inner wall to avoid separation. It is recommended also that the 
present cae contour be cut off at the station where the value of ¢ 
is a minimum. Addition of static pressure taps on the inner wall at 
the inlet plane and a method for surveying the exit plane would pro- 
vide data to better describe the flow. 

It is recommended that the method of measuring pressure be investi- 
gated and improved to give more precise results for small, fluctuating 
pressure measurements, prior to conducting further tests on diffusers 


with Ac, approaching zero. 
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APPENDIX A 


DETERMINATION OF DIFFUSER EFFICIENCY 
The diffuser efficiency was determined by the method of Vavra [10] 
for an adiabatic compression process with friction. The polytropic 
efficiency is defined for the differential compression from p to 


p # Ap of Figure 3 as 





le arr: (Al) 


The dimensionless flow function oA and the relation for the poly- 


tropic exponent, n, from Ref. 10 are 





Ke L 
d= w T —_—— = (A2) 
s Eo l PIP > 
and 
Yq 
= A3 
Ley Cee (A3) 


where w, T, ; T)> Po and Pp, are either measured or computed for the 
particular diffuser design. 
For an incremental change in pressure from P, to py» the pressure 
P» may be expressed as p, f Ap, where Ap is small with respect to Py: 
The pressure ratios in (A2) can be expressed as 
e+e (A4) 


By the binominal expansion theorem, 
: 2/n 
(+ -+2 a hae a) epee (A5) 
Py Py nN p 
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and similarly 


7P 
(*) = 14 Bey... (A6) 
n Py 


and by assuming Ap to be small, the (Ap/p,)* and higher order terms are 
much smaller than Ap/p, and are omitted. 


From (A2) 


Cc 


5 * Eo} = “1 (22) : es ; (A7) 


Substituting expressions (A5) and (A6), and solving for n, gives 





#he tea] 
fal 
(A8) 
ee z= —t14+1+4 oP 
2 Py 
Equating (A3) and (A8) gives 
tof tt 
as T, : 1 
n 52 O2) 
— - Ap Re 
y~ ea I) = ez L) - Z 
with 
T. 
ee ae (A10) 
F aaa | 
1 
the polytropic efficiency is 
ae ve 
2 1+ (7-1) 
Py 1 
OM Paeangee (Al1) 
Ze [ Ap, | 
oe yM, 2 = =F 1 


1 
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APPENDIX B 


DESIGN OF DIFFUSER MODEL (2) 

From the results of Ref. 2 and the test data of Model (1), a design 
point of 2 = 6 and L/AR = 6 was chosen from Fig. 18 for Model (2). The 
inlet area of Model (2) was the same as for Model (1). 

From equation (42) of the main text, 

ree 7 oe = 2 ce 


or, from equation (40), and a value of Ac, = 0.1, 


ae 
E GQ) ir orts, (B1) 


For the range of oe from 0.65 to 0.72 for a L/AR = 6, equation (B1) 
gives 


ot 
“"T-@ite.) i 
pr 


giving a range of A, from 39 to 51 nee 


2 
An approximate radius was chosen to give a reasonable axial opening 
at a radial exit plane for Ay = 45 ig: 


Since the outer surface, E, of Fig. 22 was most critical from the 
standpoint of flow separation, it was designed first. A smooth curve 
of the desired shape was sketched in for surface E. To approximate 
the curve sketched, a series of circular arcs were constructed which 
provide for simple tooling layouts. 

Figure 23 shows the method by which two circular arcs are fitted to 
satisfy the condition that both have centers on the same line at their 


tangent point (4), and that they fit into a given length, L. 
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FIGURE 22 


Design of Model (7°) 
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FIGURE 23 


Fitting of circular arcs 
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The largest radius, R is chosen to give the desired effect of 


1? 
the original sketch. A curve is then drawn from point (1) to point 

(2), where (2) is an approximate point of the desired curvature change. 
Point (2) defines a and b for any given distance d. The center of Ro 


becomes point (3). 


From geometry, R, and R‘, can be defined as 


2 2 


R, = oo (B3) 
: I - (a/R,)° 


and 


OE iy 5 5 (B4) 


Point (2) may be iterated until Ry = R',. The centers of Ry and Ry 
will then be on the same radius at their tangent point and the curve 
will fall within the axial length L and height d. 

The Fg acme as described above was extended to surfaces of 
more than two circular arc segments, such as surface E, of Fig. 19, by 
locating the are centers E-l, E-2, and E-3. Surface F was constructed 
in a similar manner. 

With areas Ay and Ay known, a linear distribution of flow area from 
station (1) to station (2) was chosen. Circles of the diameter which 
gave the desired area at a given meridional station were then con- 
structed along an approximate mean streamline as shown by the example 
circles on Fig. 22. Surface B was then sketched in, tangent to the 
area circles, and fitted with circular arcs as described above. The 


flow channel was then checked for uniformly increasing area as a function 


of length, and adjusted to prevent any convergence in the flow area. 
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APPENDIX C 


COMPUTER PROGRAM OMEGA 
Program "Omega" is a modification of program "Diffuser" used by 
Beck [6]. The program is written in Fortran IV for the IBM 360-67 
digital computer. The program was established to compute the values 
of flow rate, Re, 2, 2 ey and ue by the methods described in the main 


text. Table Cl is a listing of the program. 


Cl. Assembly of Input Data 

The field specifications for the input data Bis shown in Fig. 24. 
Card one of the data deck is a 72 space "'Hollerith" statement for data 
identification. The content of this card is printed immediately below 
the output title. Card two is a two digit I-format number indicating 
the number of data points in the deck. Card three is the input of 
inlet and exit diameters and the value of ). Cardsfour and five con- 
tain the data recorded for one data point and must be repeated for 
each point. No other cards need be repeated. 

All floating point numbers may be placed anywhere in their field. 
I-format numbers must be right justified. 


The definitions of the program input variables are given below. 


Fortran 

Symbol Definition Format 
NUM Number of sets- of data points eZ 
Dol Diameter of outer wall station one F 10.0 
DIl Diameter of inner wall station one F 10.0 
DO02 Diameter of outer wall station two F 10.0 
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Fortran 


symbol 


DO2 
DI2 
OMEGA 
IRUN 
MO 
IDAY 
PATM 
PINOZ 
DPNOZ 
TNOZ 
TEMP 
Pl 

Ql 


TTl 


TRM 


C2. Output 

The program output is in three sections. 
printout of data as read from the input cards as a check. 
section is the calculation of flow rate which is divided into input 
values converted by the program to dimensionally correct units, and 


output of flow rates computed by the program. 


Definition 

Diameter of outer wall, station two 
Diameter of inner wall, station two 
Value of shape parameter 

Number of run 

Month 

Day 

Atmospheric pressure 

Static pressure upstream of nozzle 
Pressure differential across nozzle 
Total temperature upstream of nozzle 
Temperature read at barometer 
Static pressure at station one 
Dynamic head at station one 


Total temperature measured upstream 
of station one 


Ambient room temperature at exit 


Format 


F 10.0 


FVEO.0 


F 10.0 


F 10.0 


F 10.0 


Faloe¢ 


F 10.0 


F 10.0 


The first section is 


The third section is 


also divided into input values converted by the program and output 


values of interest in evaluation of diffuser performance. 


Table C2 is a sample output of the program. 
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The second 





CARD 
a 
Data Identification 





CARD 3 














30 
DOl DO2 





Di 





CARD 4 


CARD 5 


Pl onl sige TRM | Blank 





FIGURE 24 


Input field specifications for Program Omega 
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TABLE Cl 


Listing of Program Omega 


TOS aeneowe oon iNtTS NEEDED FOR COMPUTATIONS 
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APPENDIX 1) 


REDUCTION OF DATA FROM SOVRAN AND KLOMP, REFERENCE 2 


Data of Ref. 2 were reduced by a Fortran IV program to a form 


useable for this study. 


Table Dl is a listing of the program. The 


diffuser description of Ref. 2 was the same as in Fig. 10. 


The definition of the input variables are 


Fortran Symbol 


A LFAH 

ALFAT 

DIMLL 

ARATIO 
RRATIO 
RTL 


CPR 


Descriperen 


The divergence angle, Ch» of the inner wall 
The divergence angle, Cp of the outer wall 
The dimensionless length, L/AR. 

[Meee 902 Je 105) - A,/A, 

Radius ratio, RyRy 

Radius of outer wall, station one 


Coefficient of pressure recovery 


The axial length was computed from the geometry and a numerical 


integration was made to obtain Q. The ideal pressure recovery 


coefficient was computed from the area ratio and used to compute Qe 


f° 


The results are contained in Table D2, where 


Fortran Symbol 


OMEGA 
CFBAR 
DIMLL 


OCF 


Description 


The shape parameter, 


Average local coefficient of friction, Ce 


The dimensionless length, L/AR 


The product, @ 
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TABLE D2 
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TABLE D2, continued 
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TABLE D2, continued 
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